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We describe the influence of hard wall confinement and lateral dimension on the low temperature 
transport properties of long diffusive channels and ballistic crosses fabricated in an InSb/In^Ali-xSb 
heterostructure. Partially diffuse boundary scattering is found to play a crucial role in the electron 
dynamics of ballistic crosses and substantially enhance the negative bend resistance. Experimen- 
tal observations are supported by simulations using a classical billiard ball model for which good 
agreement is found when diffuse boundary scattering is included. 



PACS numbers: 73.23.Ad 
Keywords: InSb, ballistic transport 

I. INTRODUCTION 

The InSb two dimensional electron gas (2DEG) is at- 
tractive for room temperature (RT) applications such as 
high speed logic devices [1] and high spatial resolution 
magnetic field sensors [2] where carrier mobility plays an 
important role. Recent improvements in the growth of 
InSb/In^Ali-^Sb quantum wells (QWs) on GaAs sub- 
strates have lead to RT electron mobility values in ex- 
cess of fx = 6 m 2 /Vs approaching the phonon limited 
value of 7 m 2 /Vs.[3] For applications requiring high spa- 
tial resolution, device miniaturization inevitably leads to 
the relevant lateral dimensions of the conducting channel 
becoming comparable to the elastic mean free path Ao, 
where transport is ballistic and bulk properties are no 
longer preserved (kp is the Fermi wavevector). There- 
fore, it is essential to understand how the device prop- 
erties are altered when fabricated at the nanoscale. For 
example, the mobility in long InAs/AlSb 2DEG chan- 
nels fabricated using reactive ion etching (RIE) is de- 
graded from that in the bulk due to top surface damage 
caused by energetic ions, but the RIE-induced sidewall 
roughness degrades the mobility further as the width of 
the channel (w) is reduced below Ao owing to electron- 
boundary scattering. [4] Degradation of \i is detrimental 
to the performance of transistors, conventional Hall, and 
extraordinary magnetoresistor (EMR) sensors based on 
diffusive transport, but it is not clear how properties are 
further effected in the mesoscopic regime. 

When the length of the channel (I) is reduced below the 
mean free path (I <Ao), electrons can traverse the device 
without scattering internally and the channel resistance 
is expressed in terms of the transmission probabilities 
between reservoirs attached to each lead, following the 
Landauer-Bttiker (L-B) formalism. [5] Ballistic transport 
in GaAs/ALcGai-^As micro junctions (where l,w > Ao) 
has been widely studied at low temperatures and a good 
understanding of the phenomena is established. [6,7] A va- 



riety of distinct departures from classical behaviour ap- 
pear in the low field magnetotransport of simple cross 
junctions, such as a negative resistance in zero magnetic 
field referred to as bend resistance, [8,9] and a quenched or 
negative Hall resistance at low fields. [10] The above men- 
tioned anomalies can be adequately described in terms of 
classical electron trajectories by treating the electrons as 
classical particles which, in analogy to ray optics, reflect 
from the boundaries with predicable trajectories. [11] Lat- 
eral depletion of conducting channels, or sidewall deple- 
tion, is also relevant as devices are miniaturized as this 
limits the minimum device dimensions. With the excep- 
tion of the InAs system that exhibits very little side- 
wall depletion, [12] Fermi level pinning at the surface of 
mesa etched III-V devices can lead to substantial side- 
wall depletion, which is not straightforward to deduce, 
however, knowledge of the depletion width (wdep) is es- 
sential in order to determine the true effective electrical 
width (w e ff) of narrow channel devices e.g. sub-micron 
Hall sensors[13] and quasi-lD wires[14]. 

Experimental knowledge of the mesoscopic properties 
of InSb and its heterostructures is still relatively lim- 
ited. 15 Negative bend resistance (NBR) was reported 
InSb/In^Ali-^Sb sub-micron structures up to T < 205 
K.[16] It was proposed that parallel conduction in the 
heterostructure masks the ballistic component from the 
2DEG Indeed, a recent study of transport in similar 
InSb/In^Ali-^Sb samples showed that at elevated tem- 
peratures intrinsic conduction in the ternary buffer layer 
contributes up to ~ 5% of the total conduction. [17] The 
significance of such parallel conduction is accentuated in 
shallow etched sub-micron structures. This technologi- 
cal problem may be overcome by improved heterostruc- 
ture design. Therefore, two regimes are identified in 
InSb/In^Ali-^Sb sub- micron structures (a) low temper- 
atures (< 100 K), where ballistic transport in the 2DEG 
is dominant and (b) high temperatures (> 150 K) where 
as yet, in all reported structures, parasitic intrinsic con- 



duct ion in the buffer layers occur. 

We emphasise that the interaction of charge carriers 
with the device boundaries plays a central role in deter- 
mining the characteristics of sub-micron devices; in par- 
ticular, ballistic anomalies are acutely sensitive to the de- 
vice dimension, geometry, [18] and the specularity of the 
boundary scattering. [19,20] Accordingly, we report here 
a detailed study of the influence of device size, sidewall 
depletion, and boundary scattering on the magnetotrans- 
port properties of InSb/In^Ali-^Sb mesoscopic struc- 
tures with hard wall confinement. For the purpose of 
this article, we present data from long channels and sub- 
micron crosses with lateral dimensions down to w ~ 170 
nm, and we restrict ourselves to low temperatures where 
intrinsic conduction is negligible. A detailed analysis of 
the ballistic transport anomalies and the agreement with 
theory is presented with the aid of a classical billiard ball 
model. 



II. EXPERIMENTAL METHODS 

Devices were fabricated from a single modulation 
doped InSb/In^Ali-^Sb QW heterostructure grown by 
molecular beam epitaxy onto a GaAs (001) substrate. In 
growth sequence, the sample consists of an AlSb (200 
nm)/Ino.gAlo.lSb (3 fim) buffer layer, a 30 nm InSb QW 
and, a 50 nm Ino.85Alo.15Sb cap in which a single Te 
^-doping layer is located, 20 nm above the top of the 
QW. The properties of the as-grown 2DEG were deter- 
mined from a 40 fim wide Hall bridge (control sample) 
fabricated using conventional wet etching. At 2 K the 2D 
electron density (n) and mobility (/1) were n = 3.95x1015 
m -2 and /i= 19.5 m 2 /Vs, corresponding to a mean free 
path of Ao= 2.03 /am and a Fermi wavelength of 40 nm. 
The bulk magnetotransport properties of this and sim- 
ilar samples were recently reported.3,17 Measurements 
were performed with the sample in the dark using a low- 
frequency lock-in technique (currents between 100 and 
500 nA) and with magnetic field B applied perpendicu- 
lar to the plane of the 2DEG. 

Hall crosses and Hall bridges with varying w were pat- 
terned by electron beam lithography using negative tone 
resist as an etch mask. Pattern transfer was achieved 
using an inductively coupled plasma- RIE in a CH4/H2 
gas mixture at a pressure of 10 mTorr, forming shallow 
mesas of w 135 nm depth that provide hard wall con- 
finement. The process parameters yielded an etch rate 
of the ternary Ino.85Alo.15Sb compound of ~ 10 nm/min. 
Ti/Au Ohmic contacts were made using standard optical 
lithography and a cold shallow contacting technique. 21 A 
deep wet chemical etch was used to remove the entire 3 
jam thick buffer layer surrounding the device and contacts 
the volume of remaining buffer layer beyond the shallow 
boundaries of the crosses was minimised by mask design 
and controlled lateral etching [see Fig. 1(a)]. Electron 
micrographs of a w = 171 ± 10 nm cross and w = 550 ± 
10 nm Hall bridge are shown in Fig. 1(a) and (b) [the un- 
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FIG. 1. Electron micrographs of (a) a typical device structure 
showing the leads and contact arrangement for a cross, and 
(b) a 550 nm wide Hall bridge. Inset to (a): A w=171 nm 
cross, (c) A schematic of the Hall cross geometry. 



certainty in w is due to residual polymer deposit from the 
RIE at the mesa edge (fencing)]. The junction corners 
are nominally square, but a small unavoidable rounding 
of the corners results from the large proximity effect in 
the e-beam lithography of InSb. 



III. CLASSICAL BILLIARD BALL MODEL 

We calculate the bend and Hall resistance of the cross 
junction following the classical model of Beenakker and 
van Houten that treats electrons as classical particles (bil- 
liard balls) reflecting from the device boundaries. [11] The 
resistance in the ballistic regime is expressed in terms of 
the transmission probabilities between the various leads 
by the L-B formula. [5] We consider the four-terminal 
hard-wall cross geometry with four-fold symmetry, in 
which case, respectively the Hall and bend resistances, 
RH and RB are given by: 
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where Tp, and T# are the probabilities of an elec- 
tron transmitted from the injection lead (arbitrary) to 
the forward, left and right hand leads respectively, and 
Ro = /i/2e 2 7Vwith TV equal to the number of transverse 
modes at the Fermi energy. In the semiclassical limit and 
for hard wall confinement, N is given by N = kpui/it. 
All calculations presented are for N ^> 1. The geom- 
etry of the cross junction is shown in Fig. 1(c) and 
is defined by three parameters: the lead width w, lead 
length T and radius of curvature of the corners, r, with 
r 2 = x 2 +y 2 in the plane. The transmission and reflection 
coefficients are calculated by injecting a large number of 
classical particles (5xl0 4 ) from a specified injection lead 
uniformly across the lead with an angular distribution 
P((j)) = l/2cos(0) ((/> being the angle with respect to the 
lead axis). [11] The trajectories of the particles are deter- 
mined via integration of the equations of motion using the 
Verlet technique until they exit the junction via one of 
the four leads. Particles are injected into the junction re- 
gion at the Fermi velocity vp = Tikp / m* with an effective 
mass m* which takes into account the modifications due 
to band non-parabolicity within an analytical model for 
the dispersion, E(l + aE) = h 2 k 2 7 /2m* where a is non- 
parabolicity parameter. [22] For the InSb QW studied 
here we use a subband edge effective mass m* = 0.0162 
and a non-parabolicity parameter of a = 3.8eV _1 which 
gives a fit to an 8 band k.p model of a 30 nm QW with 
Ino.85Alo.15Sb barriers to within a few meV over a 100 
meV range. 

We incorporate diffuse boundary scattering into the 
model using the approach of Blaikie et al. [20] . Boundary 
scattering is captured using a single specularity param- 
eter, p, that describes the probability of a particle scat- 
tering diffusively (1-p) from a boundary. After a diffuse 
scattering event, particles are re-injected at the collision 
point with an angle —it/ 2 < < tt/2 from the bound- 
ary normal chosen randomly from a uniform distribution. 
Within this model, the transmission coefficients are sen- 
sitive to the lead length V as this directly affects the num- 
ber of interactions with the boundary. 



IV. RESULTS AND DISCUSSION 
A. Diffusive properties of long channels 

Fluctuations in the electrostatic potential profile of a 
conducting channel can alter the transport properties via 
electron-boundary scattering, particularly in sub-micron 
devices where the channel width w < Ao and electrons 
can travel ballistically between the channel boundaries. 

Electron-boundary scattering can be characterised by 
two parameters; the specularity parameter p and A#, 



TABLE I. Properties of the InSb 2DEG obtained from Hall 
bridges with varying physical width, w at 2 K. Data for w=40 
/im represents the control sample. ^ calculated using the ef- 
fective electrical width determined in Section IVB. 
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the average distance an electron travels before the prob- 
ability of it scattering diffusely is equal to one. [19] In 
general p < 1 for both mesa etched and split-gate de- 
vices. [19,23] Xb is proportional to w, such that as w is 
reduced the electron-boundary interactions manifest in 
the transport properties. The increased backscattering 
in narrow channels enhances the zero-field longitudinal 
resistance R xx (0), resulting in an effective /i that is re- 
duced from that in a wide sample. For partially diffuse 
scattering (p < 1), a distinctive low field peak appears in 
the R XX (B) (discussed in Section IVC).[7] 

Measurements were performed on long channels in the 
Hall bridge geometry (Fig. 1(b) and inset to Fig. 2) 
with a longitudinal voltage lead separation of I = 8.4//m 
(> Ao) ensuring that transport is diffusive along the chan- 
nel. In Fig. 2 we show the longitudinal R xx and trans- 
verse R xy resistance as a function of magnetic field at 2 K 
for a 3 /im and 550 nm wide Hall bridge. Shubnikov de- 
Haas (SdH) oscillations in R xx are observed in each de- 
vice superposed onto an increasing background resistance 
related to parallel conduction in the upper barrier. [17] 
The 2D electron density n is determined from the peri- 
odicity of SdH oscillations and the mobility /i from the 
zero field resistance, according to /1 = l/wR xx (0)ne. The 
experimental n, \i and the corresponding mean free paths 
for the w = 3 /im and 550 nm Hall bridges are given in Ta- 
ble I, together with the properties of the control sample 
(w = 40 /im). A monotonic decrease in n and /i is ob- 
served as w is reduced. The reduction in n is attributed 
to the lateral potential formed by a sidewall depletion 
region (discussed further in the section IVB) in addition 
to the lateral confinement imposed on narrow channels 
which raise the conduction band edge in the centre of the 
channel as w is reduced, hence depleting the 2DEG. The 
observed degradation of /i is consistent with the pres- 
ence of boundary scattering which becomes increasingly 
important as w is reduced, as discussed. However, the 
mobility in the 550 nm wide channel is only approxi- 
mately 25 % smaller than in the control sample with a 
corresponding mean free path of Ao = 1.5/im. Ballistic 
transport is therefore expected in the sub-micron crosses 
(Aq > w) discussed in Section IVD. 



B. Determination of depletion width 

An important parameter of narrow channels is the elec- 
trical width w e ff. Due to the Fermi energy pinning in 
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FIG. 2. Longitudinal i? xcc (left axis) and transverse R xy (right 
axis) magnetoresistance of w=550 nm (solid lines) and 3 fim 
(dashed lines) Hall bridges at 2 K (/ = 8.4 /mi). Inset: A 
schematic of the device structure and the relevant dimensions. 
The positions of Bm ax and B m in relate to features associated 
with boundary scattering (see text). 



the band gap at the air-interface, sidewall depletion is fre- 
quently observed for narrow mesa-etched channels result- 
ing in a w e ff that can be substantially smaller than the 
physical width, [24,14] The difference is equal to the 
lateral depletion width at each boundary, 2wd e p- Knowl- 
edge of Wdep is essential for many applications but is 
not straightforward to gain. We found that devices with 
w < 134 nm were electrically depleted over the entire 
temperature range. This puts an initial estimate on the 
depletion width at Wd ep ~ 67 nm. We determine w e ff 
from tracking the depopulation of quasi- ID magnetoelec- 
tric subbands in the low field R xx data of narrow chan- 
nels. [25] Like 2D Landau levels, these hybrid subbands 
depopulate with increasing field, but do so at a slower 
rate, as is evidenced by a non-linear subband index (i) 
versus 1/B plot. For a parabolic confining potential, the 
magnetic depopulation of subbands is described by [25] 



— N 1d uj I — 



(2) 



where Nid is the ID electron density, ujq is the character- 
istic frequency defining the strength of the confinement 
and uj = (uJq + uj^) 2 where uj c = eB/m* is the cyclotron 
frequency. One can see that for small fields, the depen- 
dence ofi on 1/B is non-linear and for large fields, uo — >> u c 
and i is proportional to 1/B as in the usual 2D case. A 
subband depopulation diagram for the w = 550 nm Hall 
bridge is shown in Fig. 3. A pronounced departure from 
linear behaviour (dashed line) is observed below 1 T al- 
lowing us to implement the model of Ref. 30. The solid 
line in Fig. 3 represents a least squares fit to the data 
using Eq. 2 with an effective mass at the Fermi energy 
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FIG. 3. Subband depopulation diagram for the iu=550 nm 
Hall bridge at 2 K. The faint solid line shows the correspond- 
ing R xx data (right hand axis) from which the subband indices 
(left hand axis) were assigned. The dashed and solid lines rep- 
resent fits of Eq. 2 to the high field linear portion of the data 
and the low field non-linear respectively. 



of m* = 0.022mo (we found that the fitting results are 
relatively insensitive to small changes of w 10% in m*). 
The confinement energy Tiuoq and N\d are determined to 
be 2.6 meV and 3xl0 9 m _1 , respectively. The effective 
width is then estimated from[25] 
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Substituting the values of ujo and N±d into Eq. 3 we de- 
termine w e ff = 414 nm ± 5 nm. This implies a depletion 
width of w^p = (w — w e ff)/2 = 68 nm ± 6 nm which 
is in remarkably good agreement with the estimate made 
directly from the electrical depletion of devices of w < 
134 nm. 

Finally, we remark on a separate and consistent es- 
timate of w e ff made from a classical size effect. The 
electron backscattering in narrow channels that enhances 
Rxx(fyi is suppressed by a perpendicular magnetic field 
due to the formation of localised edge states, or classical 
skipping orbits at the boundaries. This leads to a nega- 
tive MR peaked at B = 0, persisting until B min = 2Bq, 
where Bo is the field when the cyclotron radius, R C1 equals 
w e ff, at which point a marked change in slope is ex- 
pected. [26] As seen in Fig. 2 (and more clearly in Fig. 
4), this behaviour is observed in our data. A kink in the 
low field MR is observed at a field B m i n ^ 0.5 T (indi- 
cated by an arrow) , from which we estimate w e f / ~ 406 
nm (i.e. Wd ep ~ 72 nm). This estimate is consistent 
with the value obtained from the magnetodepopulation 
analysis, adding confidence to our estimate of Wdep- 
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C. Partially diffuse boundary scattering in narrow 
channels 

The specularity of the boundary scattering plays a cru- 
cial role in the transport of submicron devices. In par- 
ticular, Blaikie et al [20] showed that resistance anoma- 
lies in ballistic devices can be substantially enhanced by 
partially diffuse boundary scattering. The specularity 
of boundary scattering can be studied from measure- 
ments on long narrow channels (w < Aq) where electron- 
boundary interactions manifest in the resistance. It has 
been shown that partially diffusive boundary scatter- 
ing leads to an anomalous peak in R xx at small fields 
(0 < B < B m i n ) with a position (B max ) that scales in- 
versely with u>. [19,27,28] As seen in Fig. 2 a pronounced 
peak is distinguished in the low field R xx of the w = 550 
nm Hall bridge at B max « 180 mT (indicated by the 
arrow). We note that a low field peak with entirely dif- 
ferent origin was also predicted [11] and experimentally 
observed [29] in the MR of ballistic Hall bridges where 
/ < Ao- In our case,/ ^> Ao so that the measurement is in 
the diffusive regime and the observed peak is unambigu- 
ously attributed to partially diffuse boundary scattering. 

In Fig. 4 we show the low field MR of the w = 550 nm 
Hall bridge [R xx (B)—R xx (0)]/R xx (0) plotted against the 
normalised field B /B (using w e ff = 414 nm) at various 
temperatures between 2 K and 80 K after subtraction 
of the high field quasi-linear background. The classi- 
cal model for in-plane MR of thin metal films (where 
the film thickness t < Ao and p = 0) predicts that 
Bmax = 0.55B i.e. when w e ff/R c = 0.55. This has 
been considered as a method of estimating w e ff. [27-29] 
We have found that the boundary scattering peak oc- 
curs at a somewhat larger value B max « 0.7Bo- It fol- 
lows that estimating weff from the classical prediction [27] 
w e ff = 0.55R C yields a value significantly less than that 
obtained in the previous section. Given that the calcula- 
tions of B max are sensitive to the details of the model [29] 
and that predicted values have been reported in the range 
0.55Bo < B max <Bo, [29,30] we suggest that this method 
provides a less reliable estimate of w e ff. 

The decay of the peak at B max with temperature is 
associated with the reduction of Ao in the bulk of the 
channel below A# (taken from the control device where 
boundary scattering can be neglected). [19] The bound- 
ary scattering length is estimated (rather arbitrarily) by 
assuming that A# ~ Ao (in the bulk of the channel) at 
the temperature (T f ) when AR xx (B rnax )/R xx (0) = l.[23] 
Taking V « 40 K corresponds to \ B ~ 1.75/im. The 
specularity parameter p is then estimated from the em- 
pirical relationship 1 — p « w e ff/\B yielding p w 0.71 
for the w = 550 nm (w e ff = 414 nm) Hall bridge. The 
inset of Fig. 4 shows the amplitude of the peak at B max 
plotted against Ao obtained from the control sample at 
each temperature. The amplitude was extracted with 
respect to a straight line drawn between data at B /Bo 
= and B /Bo = 2. Using this plot, A# may be inter- 
preted as the value of Aq when the peak amplitude de- 
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FIG. 4. The magnetoresistance, AR xx (B)/R xx (0), of the w 
= 550 nm Hall bridge plotted against the normalised field 
B/Bo, at various temperatures after subtraction of a linear 
background. Inset: The peak amplitude plotted against the 
mean free path (Ao) in the control sample at each tempera- 
ture. 



cays to zero. Two dependences on Ao are distinguished 
in the data, a rapid decay (solid line) labelled as A and 
a slower decay (dashed line) labelled as B. We broadly 
separate these into the regimes where remote ionized im- 
purities and phonons dominate momentum scattering in 
the bulk of the channel, respectively. We consider regime 
B unsuitable for this analysis since large angle phonon 
scattering randomises the electrons momentum in addi- 
tion to diffuse boundary scattering events which alter A#. 
Therefore, only at low temperatures (regime A) can infor- 
mation on Ab be extracted with confidence. In regime A, 
we extrapolate a value of A# ~ 1.85/im, giving p « 0.77 
which is similar to the previous estimate. We conclude 
from our analysis that p ~ 0.7 — 0.8. 

The value of p is expected to be a property of the 
boundaries themselves and therefore be the same for de- 
vices fabricated in the same way. Given the assumptions 
made to estimate p, emphasis should not be on the value 
of p itself but rather it should be sufficient that one ob- 
serves the characteristic low field MR features shown in 
Figs. 2 and 4, to conclude that partially diffuse boundary 
scattering is significant and p < 1. 



D. Ballistic transport in cross junctions 

We now turn to the experimental results in ballistic 
crosses formed from two intersecting channels of width 
w [see inset to Fig. 1(a)], where the relevant lateral di- 
mensions L « w are substantially less than the mean free 
path. We present the results from four crosses with phys- 
ical widths (inferred from SEM inspection) of w = 924, 
550, 400, and 171 nm ± 10 nm. The inferred effective 
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FIG. 5. The Hall resistance Rh in units oih/e 2 as a function 
of B for a w = 171 nm (blue line), 400 nm (black line), 550 
nm (red line), and 924 nm (green line) cross. The dashed 
lines represent the classical 2D result. Top inset: Low field 
data illustrating the anomalies in Rh- Data for the w — 
550 nm and 924 nm crosses are offset by 0.5 kQ for clarity. 
Bottom inset: Dependence of n on the inferred effective width 
w e ff = w — 2wde P of the devices (crosses and bridges). 

electrical widths w e ff = w — 2wd e p are given in Table 
II where we have used the depletion width determined 
in Section IVB (w^ep = 68 nm). Note that the smallest 
cross (w — 171 nm) has an estimated electrical width of 
w e f f w 35 nm which is among the narrowest conducting 
mesa-etched devices reported. [12] 

Figure 5 shows the results for the Hall resistance 
Rh = V^2/h,3 [the lead arrangement is shown in Fig. 
1(a)] as a function of B for crosses. Here Vij and J mn 
indicate the voltage of terminal i measured with respect 
to j when current is passed from terminal m to n, respec- 
tively. Data for w = 924 nm, 550 nm, and 400 nm were 
taken at 2 K and the w = 171 nm at 40 K (the w — 171 nm 
junction became depleted for T < 30 K). Quantum Hall 
plateaus are resolved in the data from the largest three 
crosses. The electron densities are determined from fits 
of the classical 2D result Rh(B) = —B/ne (indicated 
by the dashed lines in Fig. 5) to the high field linear 
portions of data. Rh(B) for the w = 171 nm cross is 
strikingly different - no obvious quantisation of Rh oc- 
curs over the entire field range and Rh(B) is non-linear 
up to | B | ^4T making the determination of n less triv- 
ial. The extracted electron densities of the crosses are 
listed in Table II. The dependence of n on w e ff is pre- 
sented in the bottom inset to Fig. 5 including data from 
the wider Hall bridges. 



At low fields \B\ < 1 T, clear anomalies appear in Rh 
for all crosses (top inset to Fig. 5) - the development 
of the anomalies with decreasing w is clear. No sup- 
pression of Rh around B = is observed in the largest 
three crosses (a very small reduction is found for the w 
= 400 nm cross). In the smallest cross (blue line) the 
effect is striking; Rh is completely quenched and nega- 
tive (positive in our configuration) up to \B\ < 0.65 T. A 
small asymmetry in Rh(B) is observed in all cases which 
is attributed to asymmetries in the geometric junction. 
The appearance of quenching is of interest with respect 
to the geometry of the junction. Baranger and Stone [31] 
showed that generic quenching of Rh only occurs in junc- 
tions with rounded corners; a consequence of a horn colli- 
mation effect [32] which yields a non equilibrium momen- 
tum distribution that enhances the forward transmission 
(Tp) at the expense of the transmission into the left (Tc) 
and right (Tr) leads [c.f. Eq. 1(a)]. Electron collima- 
tion was experimentally verified by Molenkamp et al. [33] 
and is a key concept in describing ballistic anomalies as 
we demonstrate here. Likewise, the negative Rh results 
L |rj3m rebound trajectories (directing electrons into the 
wrong lead for a given field direction) that are only ef- 
fective in rounded junctions when the radius of curva- 
ture of the junction corner (r) is large compared to the 
lead width i.e. r/w > 1. [34,18] The appearance of these 
features is therefore a clear signature of both significant 
rounding and collimation in the w = 171 nm cross. The 
former is perhaps surprising because the junctions are 
nominally square, however, we accept some small round- 
ing is inevitable in the e-beam and etch process which 
emphasises that w e ff must be very small in this case. 
Conversely, the lack of quenching in the largest three 
crosses implies relatively little collimation and small r/w 
i.e. the junctions are approximately square. 

Beyond the quenched region, Rh rises above its classi- 
cal value (dashed lines in Fig. 5) in all devices, marking 
the onset of the classical last plateau. [10,18] At larger 
fields still, the data rejoins the classical Hall resistance 
(indicated by the arrows in Fig. 5). For the w = 171 
nm cross the non-linearity persists up to \B\ ~4T. The 
sharp rise in Rh(B) above its classical value results from 
trajectories that guide electrons into a side lead with min- 
imal boundary reflections thereby enhancing the asym- 
metry between Tl and Tr.[11] When guiding is fully ef- 
fective, electrons are no longer reflected back into the 
junction (skipping orbits along the junction perimeter) 
and T F ,T R <C T L « 1. With reference to Eq. 1(a), 
in this regime Rh(B) is predicted to plateau at a value 
equal to the contact resistance of the lead Ro = h/2e 2 N. 
For B > 2Bq classical behaviour is recovered. Although 
a clear plateau region is not observed in our experimen- 
tal data, features consistent with the predictions of the 
classical model are observed. For example, the estimated 
value 2B = 4.5 T for the w = 171 nm cross coincides 
approximately with the field at which the experimental 
data rejoin the classical Hall slope. Similar agreement 
is found for each cross indicating that the estimations of 
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FIG. 6. (a) Bend resistance Rb = V4, 3/^1,2 as a function of 
magnetic field for a w = 171 nm (blue line), 400 nm (black 
line), 550 nm (red line), and 924 nm (green line) cross. The 
w = 171 nm data is plotted on a difference scale for ease of 
comparison. The dotted curve is the reciprocal measurement 
Rb = Vi ,2/^4,3 for the w — 400 nm cross illustrating that 
asymmetries in B originate from junction asymmetry, (b) 
Dependence of the experimental NBR amplitude ARb (O) 
on 1/N = 7r/kFW e ff and w e ff. The results from billiard 
model simulations using the parameters given in Table II ( x ) 
and results for p = 1(+) are also shown. The dashed line is 
a guide to the eye illustrating a 1/N dependence. Inset: A 
schematic showing the definition of ARb- 

w e ff and kp = (27rn)^ are close to the true values (Rq 
and Bo for each cross are listed in Table II). 

In Fig. 6(a) we show low field results obtained in 
the bend resistance configuration Rb = Va,z/Ii,2 for the 
crosses (again, data for the w — 171 nm cross was ob- 
tained at 40 K). A substantial NBR peak centred about 
B = is observed in all devices that increases as w is 
reduced. Asymmetries in the field dependence are also 
observed in this configuration and are particularly evi- 
dent in the data for w = 400 nm cross [solid black line 
in Fig. 6(a)]. To ascertain the origin of the asymme- 



tries, measurements were repeated with the current and 
voltage leads interchanged. Representative data for the 
w = 400 nm cross is shown by the dashed black line 
in Fig. 6(a). One can see that the reciprocity relation 
Rmn,ij(B) = Rij, mn (—B) is obeyed demonstrating that 
the field asymmetries indeed originate from asymmetries 
in the junction geometry. [5] This is representative of each 
device measured. 

The origin of NBR is well established: it arises from 
straight through trajectories which raise the potential at 
lead 3 with respect to lead 4 [see Fig. 1(a)], resulting in 
a negative resistance. This corresponds to Tp ^> Tl, Tr 
in the L-B formula [c.f. Eq. 1(b)]. In a small magnetic 
field the Lorentz force curves the trajectories into the 
correct lead 4 and the NBR decays to zero producing a 
characteristic negative peak (for B > this corresponds 
to Tr = Tp = 0). In our case, a small diffuse background 
resistance is present ranging from 20-30 Q. Before the 
background resistance is recovered, a small overshoot of 
positive resistance is observed in each cross [indicated by 
the arrows in Fig. 6(a)] due to rebound trajectories, in 
rounded junctions, that briefly increase the transmission 
into the opposite lead. This coincides with the rise in Rh 
to the last plateau. 

The case of NBR in zero magnetic field is useful 
because the solutions to the L-B formulae are simpli- 
fied, allowing information on the transmission proba- 
bilities and collimation to be extracted. [33] At B = 0, 
T L = T R = T s and Eq. 1(b) reduces to R B (0)/R = 
(1 - T F /T S )/[4(T S + T F )]. For symmetric hard walled 
junctions with fixed geometry (i.e. fixed r/w), the 
classical model [32] predicts a universal scaling of resis- 
tance curves when normalised by Ro and Bo- In other 
words, the transmission coefficients and hence the col- 
limation are approximately the same for junctions with 
fixed geometry.^ (0) therefore scales inversely with the 
number of channels N = kpw e ff/i\ (see Table II). For 
analysis purposes, we define a NBR amplitude ARb = 
R' b (0) — Rb(0) as the difference between the interpolated 
background resistance at B = and i?#(0) (i?#(0) = 
in the billiard model) [see inset to Fig. 6(b)]. Figure 
6(b) shows the variation of ARb with 1/N on a log- log 
plot. ARb scales approximately, but not exactly, with 
1/N (compare to the dashed line). Universal scaling pre- 
dicts that the normalised resistance Rb (0)/Ro is inde- 
pendent of w e ff and kp- Accordingly, in Fig. 7(a) we 
show ARb/Ro plotted against w e ff in our devices. Re- 
markably we find that ARb/Ro is almost identical for the 
largest two crosses and thus exhibit the traits of univer- 
sal scaling. The geometries in these crosses must there- 
fore be equivalent, which is consistent with the assertion 
that these junctions are approximately square (i.e.r/w is 
small). Scaling for these data is also preserved for B > 
when Rb/Ro is plotted against B /Bo (not shown). This 
is not true for the two smaller crosses as evidenced by 
a monotonic increase of ARb/Ro with decreasing w e ff, 
indicating that collimation in the crosses is increased as 
weff becomes smaller. These observations provide valu- 
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able insight into the geometry of the junctions which is 
used in the billiard model calculations presented in sec- 
tion IVE. 



E. Simulation results 

To explore further the electron dynamics within the 
cross junctions, calculations of the bend resistance were 
performed using the classical model described in Sec- 
tion III. Classical and quantum mechanical calculations 
of ballistic anomalies in micro junction have previously 
been performed by various authors [6, 11, 29, 35] and as 
discussed in these works, the geometry of the junction 
determines the magnitude and character of the resis- 
tance anomalies. The parameters in the calculations are 
w,r,l f ,vp, and p (a schematic of the cross geometry is 
repeated in Fig. 7(a) inset for clarity), vp is set by the 
experimentally determined kp (see Section III). We start 
by considering the magnitude of the experimental NBR 
and its implications on the collimation in the crosses, and 
then compare our results for Rb(B) with experimental 
data. 

Simulations of A.Rb(0)/Ro f° r a square (r/w = 0) 
and rounded junction (r/w = 2) with specular bound- 
ary scattering (p = 1) are shown by the solid lines in Fig. 
7(a). Recall that no collimation occurs in the square 
junction, when p = 1, whereas collimation is induced in 
the rounded junction via the horn effect. The experimen- 
tal ARb(0)/Ro of the two largest crosses (which we ex- 
pect to be approximately square) exceeds the calculated 
values for a square junction by factor of « 4 and even a 
rounded junction by a factor of « 2. The anomalously 
large NBR implies additional collimation is present other 
than the horn effect which we attribute to the diffuse 
collimation effect [20] that results from partially diffuse 
boundary scattering (p < 1) in our devices (as shown in 
Section IVC). The origin of diffuse collimation is the in- 
creased backscattering of electrons that enter the leads 
with large angles <\> with respect the lead axis. There- 
fore, electrons injected with a |cos(0) distribution arrive 
at the junction region after traversing a lead of length V 
with a distribution more strongly peaked in the forward 
direction (hence increasing the ratio Tp/Tg). The result- 
ing angular distribution differs from the horn effect result 
in that it is more sharply peaked in the forward direc- 
tion. [29] Consequently, the NBR for p < 1 has a distinc- 
tively sharper and more triangular shape about B = 
than in the p = 1 case. In support of this conjecture, the 
experimental data in Fig. 6(a) exhibit the characteristic 
sharp NBR associated with diffuse collimation. 

Billiard simulations with p < 1 were implemented us- 
ing the approach of Ref. 20 (for details see Section III). 
To illustrate the enhancement of the NBR from diffuse 
collimation a calculation of ARb(0) /Rq for a square junc- 
tion with p = 0.7 and I' /w = 3 is shown in Fig. 7(a) by 
the dashed line. Remarkably, even for a square junction, 
the NBR amplitude is increased by a factor of ~ 5 over 
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FIG. 7. (a) Normalised NBR amplitude AR B /Ro of the 
crosses plotted against w e ff. Horizontal lines represent the 
results from the billiard model for square and rounded junc- 
tions with p = 1 (solid lines) along with a square junction 
with p = 0.7 and I' /w — 3 (dashed line), (b) - (d) Com- 
parisons between experimental Rb (black lines) and billiard 
model simulations of Rb (solid red lines) for three crosses 
(N ^> 1) using parameters listed in Table II. Simulations 
with p = 1 are shown for comparison (dashed red lines). In- 
set: A schematic of the geometry used in the simulations. 
Ro = (h/2e 2 )(7v/k F w eff ). 



the p = 1 case using reasonable parameters. Diffuse col- 
limation is sensitive to the ratio V/w since this directly 
influences the number of boundary collisions. In our de- 
vices we define V as the length from the junction to the 
point at which the lead width flares out [e.g. see Fig. 
1(a)]. These values are listed in Table II. Therefore we 
can simulate the whole Rb(B) curve using experimen- 
tally determined parameters n,w e ff, and V with only p 
and r as variables (w e ff = w in the billiard model). We 
note that the position and magnitude of the overshoot 
in Rb is determined by r. This was used as a method 
of optimising r from which we found that r = 100 nm 
yielded reasonable agreement with the experimental data 
for each cross. This is consistent with the assertion that 
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TABLE II. Relevant parameters for the ballistic crosses and parameters used in the billiard calculations. Effective widths 
w e ff = w — 2wde P were calculated using Wde P = 68 nm determined in Section IVB. 



w (nm) 


w eff (nm) 


n(10 lb m- 2 ) 


N 


Ro(fi) 


Bo (T) 


r (nm) 


I' (p m) 


P 


924 


788 


3.85 


39 


332 


0.13 


100 


2.5 


0.79 


550 


414 


3.77 


20.3 


638 


0.25 


100 


1.5 


0.8 


400 


264 


3.1 


11.7 


1104 


0.35 


100 


1.2 


0.69 


171 


35 


2.25 


1.3 


9775 


2.24 




0.8 





r/w is small in the largest two crosses (from the observed 
scaling) and the fact that the unintentional rounding re- 
sults from the fabrication process that should approxi- 
mately independent of w. The value of p was used as the 
only fitting parameter to adjust Rb(0) to match the ex- 
perimental Rb(0). These calculations were performed for 
each of the three largest crosses where the classical model 
is applicable (N ^> 1). The results of the calculations are 
shown in Figs. 7(b), (c) and (d) indicated by the solid 
red lines and compared to the experimental data (solid 
black lines). The agreement with the experimental data 
is excellent considering the few adjustable parameters in- 
volved, indicating that the experimental w e ff and n are a 
close representation of the true values. The diffuse back- 
ground resistance observed in the experimental data is 
likely the result of finite momentum scattering times (r) 
within the crosses, implying that not all of the electrons 
are fully ballistic as they are treated in the model (r = oo 
in the current model). Given the agreement with the bal- 
listic model, we speculate that the momentum scattering 
in the crosses should not affect the electron trajectories 
considerably and therefore the extracted p parameters 
are meaningful. This picture is supported by recent work 
showing that at low temperatures the mobility in these 
InSb QWs is dominated by small angle remote ionised 
impurity scattering. [17,36] The values of p used for the 
three crosses lie in the range 0.69-0.8 which is consistent 
with the value p ~ 0.7 - 0.8 estimated from the measure- 
ments on long narrow channels described in Section IVC. 
This result supports the assertion that p is a property of 
the boundary and thus relatively independent of w for 
devices fabricated under the same conditions. The incor- 
poration of remote ionised impurity scattering into the 
billiard model is the subject of further work. Simulations 
for p = 1 are shown by the dashed red lines in Fig. 7(b) 
- (d) for comparison. These results illustrate the more 
rounded profile of the NBR in the p = 1 case and, more- 
over, that the incorporation of realistic partially diffuse 
boundary scattering is crucial for the accurate modelling 
of device characteristics. 



V. CONCLUSIONS 

In summary, we have investigated the varia- 
tion of the low temperature transport properties in 
InSb/In^Ali-^Sb mesa-etched mesoscopic devices with 
hard wall confinement when the lateral dimensions are 
reduced below the mean free path. Measurements on 



long channels and Hall crosses fabricated from the same 
sample show that the lateral depletion width is approx- 
imately 68 nm and that boundary scattering from the 
sidewall is partially diffuse. A specularity parameter p « 
0.7 - 0.8 was deduced. Ballistic crosses show characteris- 
tic resistance anomalies in good agreement with the pre- 
dictions of the classical model and in all cases exhibit 
a significantly enhanced negative bend resistance due to 
partially diffuse boundary scattering from the sidewalls. 
Our observations are supported by classical simulations 
of the electron trajectories in ballistic crosses which quan- 
titatively accounts for both the magnitude and width of 
the negative bend resistance, using experimentally deter- 
mined parameters, and a specularity parameter p in the 
range 0.69 - 0.8. 
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